Fano-type interference in quantum dots coupled between metallic and 

superconducting leads 



O 
(N 






i 

c3 



i 

c 
o 
o 



> 

in 

O 



% 




J. Barahski and T. Domanski 
Institute of Physics, M. Curie- Sklodowska University, 20-031 Lublin, Poland 

(Dated: August 1, 2011) 

We analyze the quantum interference effects appearing in the charge current through the double 
quantum dots coupled in T-shape configuration to an isotropic superconductor and metallic lead. 
Owing to proximity effect the quantum dots inherit a pairing which has the profound influence on 
nonequilibrium charge transport, especially in the subgap regime \V\ < A/|e|. We discuss under 
what conditions the Fano-type lineshapes might appear in such Andreev conductance and consider 
a possible interplay with the strong correlation effects. 



I. INTRODUCTION 

Heterostructures with nanoobjects (such as quantum 
dots, nanowires, molecules, etc) hybridized to one con- 
ducting and another superconducting electrode seem to 
be promising testing fields where the strong electron 
correlations (responsible e.g. for Coulomb blockade and 
Kondo physics [lj]) can be confronted with the super- 
conducting order [2(. Coulomb repulsion between elec- 
trons in the solid state physics is known to suppress the 
local (s-wave) pairing and, through the spin exchange 
mechanism, eventually promotes the intersite (<i-wave) 
superconductivity [3|. Mutual relation between such re- 
pulsion and the local pairing is however rather difficult 
for studying, both on theoretical grounds and experimen- 
tally. In nanoscopic heterostructures some of these limi- 
tations can be overcome by a suitable adjustment of the 
hybridization and the gate- voltage positioning of energy 
levels involved in the charge transfer [4]. They enable a 
controllable changeover between the Kondo regime and 
opposite case dominated by the induced on-dot pairing. 

Quantum dot (QD) coupled with the strength Tjv to 
metallic conductor (N) and with Ts to superconducting 
electrode (S) can exhibit the features characteristic both 
for the on-dot pairing and the Kondo effect (including 
their coexistence) [5| . Their efficiency depends on the ra- 
tio Ts/rAT. In the limit T$ 3> Fjv the underlying physics 
is controlled by on-dot pairing and manifests itself e.g. 
by the particle-hole splitting of the quasiparticle levels. 
On the other hand for Ts <C Tjv the strong correlations 
take over. Non-trivial aspects related to such interplay 
between the Coulomb interactions and the proximity in- 
duced on-dot pairing has been addressed theoretically us- 
ing various methods like: the mean field slave boson ap- 
proach [6(, the noncrossing approximation [7], perturba- 
tive scheme @, Q , constrained slave boson technique 
numerical renormalization group [1CH12| and other 
Ha ] . Also the cotunneling regime of a Coulomb blockaded 
quantum dot sandwiched between a normal and super- 
conducting lead, where charge fluctuations are strongly 
suppressed, has been discussed emphasizing the role of 
in-gap resonances [17J |. 

As far as the experimental situation is concerned it 
has been less intensively explored. The earliest transport 
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FIG. 1: (color online) Scheme of the double quantum dot 
in T-shape configuration coupled to the conducting (N) and 
superconducting (S) leads, where interference effects originate 
from the interdot hopping t. 



measurements for N-QD-S interface have been obtained 
using the multi-wall carbon nanotube deposited between 
Au and Al electrodes [lg]. Those investigations concen- 
trated however on the specific regime UbTk > A, when 
the Coulomb correlations dominated over the proximity 
effect. Other studies of the same group have been done 
for similar structures replacing a metallic electrode by a 
ferromagnet [l{j]. Several recent efforts focused on the 
multitcrminal structures involving two normal and one 
superconducting electrodes as useful schemes for realiza- 
tion of: the crossed Andreev reflections tunable via gate 
voltages [20|. the Cooper pair splitters [21] . and the QD 
spin valves [22] . 

Very useful understanding of a subtle interplay be- 
tween the correlations and the induced on-dot pairing has 
been gained from recent measurements by R.S. Deacon 
et al [5j. The authors have explored the subgap trans- 
port originating from the Andreev-type scattering pro- 
cesses for several representative ratios T^/r^v using the 
self- assembled InAs quantum dots deposited between the 
golden (N) and aluminium (S) electrodes. Their measure- 
ments provided the unambiguous experimental evidence 
for: a) particle-hole splitting of the subgap conductance 
of the Andreev states when Ts > Tn, and b) enhance- 
ment of the zero-bias Andreev conductance due to forma- 
tion of the Kondo resonance at the Fermi level of metallic 



lead, as has been qualitatively suggested by our studies 
jlq and also indicated by other groups |23{ . 

The present work extends our former studies by tak- 
ing into account interference effects arising from the ad- 
ditional degrees of freedom. As the simplest prototype 
for Fano-type interference [24j we consider the setup (see 
figure QJ with a side-attached quantum dot contributing 
an extra pathway for electrons transmitted between the 
metallic and superconducting leads. Our analysis is com- 
plementary to the previous study by Y. Tanaka et al who 
considered the double quantum dots coupled between N 
and S leads in a T-shape setup but assuming U\ = 0, 
U2 7^ [25| and in a series configuration [26J . 

In section 2 we introduce the microscopic model and 
briefly outline basic notes on the nonequilibrium subgap 
transport. In the next section 3 we discuss a unique 
way in which the Fano-type lineshapes might be observed 
in Andreev conductance, focusing on the uncorrelated 
quantum dots. In the last part (section 4) we discuss the 
influence of correlations at the interfacial quantum dot 
which seem to have remarkable influence on the low bias 
transport. We end with the summary and suggestions 
for the future studies. 
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II. THE MODEL 

For description of the hetero junction illustrated in fig- 
ure [T] we use the Hamiltonian 

H = Hn + Hn-dqd + Hdqd + Hs-dqd + Hs (1) 
where the double quantum dot (DQD) is represented by 

HpQp=y^ j eid^djq + U\ h lt n n + (t d\ a d 2cr + h.c.). (2) 

The energies of each [i = 1,2) quantum dot electrons 
are denoted by S{ and t stands for the usual intcrdot 
hopping. We restrict considerations of the correlation 
effects (section 4) to the Coulomb repulsion U\ between 
opposite spin electrons a =t, i at the interfacial quantum 
dot. 

The external reservoirs N and S of charge carriers 
are described by H N = J2k, a ^NC'^f^c^aN and corre- 
spondingly H s = J2k. a £kscLsCk<TS - Ek A 4r5C-k;s + 
A*c_kj,sCkts assuming the isotropic energy gap A. As 
usually £k/3 = £k/3 — Up denote the electron energies mea- 
sured from the individual chemical potentials /ip which 
become detuned [In — ^s = e V if a bias V is applied 
across the junction inducing the nonequilibrium charge 
flow I(V). Fano-type quantum interference effects origi- 
nating from the hopping t to side-coupled quantum dot 
i = 2 are discussed here assuming that only the interfacial 
quantum dot % = 1 is directly coupled to external leads 



FIG. 2: (color online) Density of states p(u)) of the interfacial 
i = l QD obtained in equilibrium conditions for ei = 0, £2 = 
0.3rjv, Ts = 5rjv and a few values of the interdot hopping t. 



In the wide-band limit approximation it is convenient 
to introduce the structureless coupling constants r^ = 
2ttJ2 l^k/3| 2( 5 (w— £k) which shall be used here as the en- 
ergy units. 

Interplay between the proximity induced on-dot pair- 
ing, the correlations and the quantum interference effects 
can be in practice detected by measuring the differential 
conductance dI(V)/dV. Particularly valuable for this 
purpose is the low voltage (subgap) regime \eV\ <C A. 
Under such conditions the charge current is provided 
by the anomalous Andreev scattering in which electrons 
from the metallic lead are converted into the Cooper pairs 
in superconductor with a simultaneous reflection of the 
electron holes back to the normal lead. On a formal level 
the resulting Andreev current can be expressed by the 
Landauer-type formula [9|, Il3| 

T a(V) = ^ f duT A (u)[f(u-eV,T)- f(u+eV,T)],{4) 



where f(u>,T) is the Fermi distribution function and the 
transmittance Ta(lS) = T 2 N |G 12 (u;)| depends on the off- 
diagonal part (in the Nambu notation) of the retarded 
Green's function (JS|) of the interfacial quantum dot. 



H/3-DQD = 2J ( V ^P divCkvP + n -' 
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III. FANO RESONANCES 

Fano resonances appear in many physical systems due 
to the quantum interference of the waves transmitted res- 





FIG. 4: (color online) The differential Andreev conductance 
Ga(V) versus the bias V revealing the quasiparicle peaks 

(near ±^ e\-\-(T s /ty 2 ) and Fano-type lineshapes (near ±£2) 
for the set of parameters used in figure [2] and i = 0.1 Tat. 



FIG. 3: (color online) Changeover of the interfacial quantum 
dot spectrum from the Fano (resonance and antiresonance) 
lineshapes to the effective four-peak structure upon increasing 
the interdot hopping t for the same parameters as in figure [2] 



simplifies to (see the appendix) 
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onantly via some discrete energy level combined with 
transmittance contributed from a continuum of other 
states. In nanoscale physics such resonances are feasi- 
ble in a variety of constructions [24|. Fano lineshapes 
are present for instance in the electron transport when 
two external electrodes are in parallel coupled through a 
quantum dot and directly via a shortcut bridge [27[ . An- 
other simple possibility takes place in the electron tunnel- 
ing using two quantum dots with considerably different 
linebroadenings |28|, [29( . In the latter case the narrower 
level is responsible for forming the Fano resonance on a 
background of the broader level. 

In this work we want to analyze similar interference ef- 
fects appearing in the anomalous Andreev current, which 
is very specific because of the particle and hole degrees of 
freedom mixed with one and other. To have a clear pic- 
ture of the underlying physics let us start by considering 
the noninteracting case Ui = when the Green's functions 
of each quantum dot can be determined exactly. 

Electron transport of the setup shown in figure [1] 
is determined by effective properties of the interfacial 
quantum dot. For this purpose we compute the ma- 
trix Green's function G(ri,r 2 ) = -T T (*(ri)* t (r 2 )) in- 



troducing the standard spinor notation &' = (dL., dij. 



and if! = (ffiy ■ In the equilibrium conditions /ijy = fig 
this function depends only on the time difference and its 
Fourier transform obeys the following Dyson equation 
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where the term X (u>) corresponds to the selfenergy of 
noninteracting case (U\=0) and S (w) accounts for the 
correlation effects (discussed in section IV). Focusing on 
the deep subgap regime |w| € Awe obtain that X (w) 



When the interference effects caused by the hopping t 
to the side-coupled QD are neglected the expression © 
becomes static (w-independent) and nontrivial physics 
of this, so called atomic superconducting limit, has been 
explored in detail by several groups @, [13, 1201 including 
ourselves [la ]. 

Taking into account the quantum interference 1 7^ we 
show in figure [2] the proximity induced on-dot pairing 
[formally arising from the off-diagonal parts of ©] illus- 
trating the energy spectrum p(u>) = --ImGn(w + i0 + ) 
obtained for strong coupling to the superconducting 
lead Ts — 5Tjy. Such coupling Ts is responsible for 
the particle-hole spli tting of the effective quasiparticle 
states formed at ±^Je\-\-(Ts/'2) 2 whereas the coupling 
Tn controls their broadening. In the particular case 
£i=0 the quasiparticle peaks appearing at ±E\ (where 
E\ = \J s'f + T's I 4) are symmetric, but for arbitrary e\ 
they are weighted by the corresponding BCS coefficients 
u 2 , v 2 = \ (1 ± e\/ Ex) [Tg|. On top of such behavior we 
clearly notice that hopping to the side-coupled quantum 
dot induces additional features appearing in the effective 
spectrum near ±£2 as the Fano resonance and antireso- 
nance. For the case of both metallic leads there would 
survive just the single Fano structure at £2 which in very 
pedagogical way has been discussed by R. Zitko [28| . 

Fano-type lineshapes (see the lower panel in figure 
are present only in the weak hopping regime t <C Tn- 
For increasing t the Fano structures gradually evolve into 
separate quasiparticle peaks illustrated in figure Phys- 
ically this can be assigned to the induced pairing on the 
side-attached QD (cfe^cfet) 7^ transmitted there indi- 
rectly via the interfacial quantum dot. Such effect again 
qualitatively differs from the structures of the DQD cou- 
pled to both metallic leads [28|, [29| . 

Interrelation between the interference and proximity 
effect can be practically investigated by measuring the 




FIG. 5: (color online) Differential conductance Ga(V) of the 
subgap Andreev current versus the source-drain bias V in a 
vicinity of the Fano structure appearing at V = ±£2/6. 



tunneling current. In figure 0] we show bias voltage 
V dependence of the differential Andreev conductance 
Ga{V) — dlA{V)/dV determined at zero temperature 
from (j4|) over a broad regime covering both the sub- 
gap quasiparticle peaks. Figure [5] illustrates the resulting 
Fano- type lineshapes Ga = G *t+i + G\ nearby — e 2 , 
where x=|eV+£2| /IV and the asymmetry parameter q 
gradually decreases upon increasing the hopping integral 
t. Our results can be thought as extension of the predic- 
tions obtained for the normal electron tunneling u sing 
the T-shape DQD coupled to both metallic leads [28l.l32l| 
onto the anomalous Andreev current where the particle 
hole mixing has the essential importance. 



IV. INTERPLAY WITH CORRELATIONS 

Coulomb repulsion between electrons of opposite spins 
can have an important influence on the spectral and 
transport properties of various nanostructures. For the 
case of quantum dots coupled to both conducting leads 
such interactions are known to be responsible for: a) the 
charging effect (if a given energy level e, is attempted to 
be occupied by more than a single electron this costs the 
system an extra energy Ui), b) the Kondo effect when 
the singlet state is formed between QD and itinerant 
electrons from the leads [3l[ . In spectroscopic proper- 
ties they are manifested by appearance of the Coulomb 
satellite around u> = Si + U and the narrow Kondo reso- 
nance at the Fermi level. For heterostructures with the 
superconducting electrodes the situation is more complex 
due to a competition between the induced on-dot pairing 
and Coulomb repulsion. 

The rich interplay between the quantum interference, 
correlations and proximity effect for the configuration 
shown in figure [T] have been so far addressed using the 
density functional technique .33] (which does not capture 
the Kondo physics) and by the numerical renormalization 
group calculations [25j ■ In latter case the authors focused 
on U\ = 0, Ui 7^ when the side-attached quantum dot 
can indirectly form the Kondo state with electrons of the 
metallic lead affecting the Andreev transport. 
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FIG. 6: (color online) Density of states p(uS) of the cor- 
related interfacial QD in the Kondo regime obtained for: 
ei = -0.75riv, T s =T N , U = 5Fiv, t = 0.1T N and temper- 
ature fcflT = O.OOirjv. The upper panel shows the spectrum 
for £2 = 0.25Fjv with two Fano structures appearing at ±£2 
aside the Kondo peak. The lower plot corresponds to £2 = 
when both the Kondo and Fano structures are superimposed. 



To account for the correlation effects predominantly 
originating from the interfacial quantum dot we extend 
here the procedure previously used by us for studying the 
single quantum dot [16j . The main idea is to approximate 
the correlation selfenergy S (w) by the diagonal matrix 
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Such assumption (applied also in the NRG studies |10| ) 
can be thought as the simplest ansatz for the many-body 
selfenergy S (ui) allowing to combine the proximity ef- 
fect © with the correlations, brought separately from 
the particle and hole channels. In more advanced treat- 
ments one should take into account the possible feedback 
effects between these normal and anomalous channels. 
We nevertheless hope that by imposing ([7]) we can get 
some insight at least on a qualitative level which might 
stimulate the future studies. 

Within qualitative studies of the correlation effects we 
can describe the Coulomb blockade and Kondo effects 
using the following equation of motion expression [31| 

w -ei - [S 3 (a;) + t/i(l-ni,oOJ 
where S I=1)3 (u;)=E k |Viciv| 2 [/(w i r)] 1 ^[(a;-Ckiv)- 1 + 



SUMMARY 
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FIG. 7: (color online) Andreev conductance Ga(V) expressed 
in units of 4e /h obtained for Ts — Tn (left h.s. panel) and 
Ts = 2rjv (panel on the right). The thin dashed lines cor- 
respond to t — whereas the thick solid lines show influence 
of the Fano-type interference for t = 0.1Fjv. In both plots we 
used £2 = and the same model parameters as in figure [6] 



(c^-C/ 1 -2e 1 4<kw)" 1 ], ni ia = {d[J la ) andw = w+-^. We 
explored the interfacial quantum dot spectrum and the 
related transport properties at fcsT = 0.001IV, i-e. well 
below the Kondo temperature. Specific numerical com- 
putations have been done for e\ — — 0.75Fjv, U\ = 5r ' n 
and symmetric coupling to both external leads giving the 
optimal conditions for any low-bias features in the An- 
dreev current [5|, [lg|. This ratio Tg/Y^ ~ 1 is a reason 
why the particle-hole splitting is hardly visible, but oth- 
erwise (for larger T s) the Kondo peak is either reduced 
or completely absent |16| . 

Upper panel in the figure [6] illustrates the Fano reso- 
nance/ antiresonance around ±£2 (where £2 = 0.25Fat) ob- 
tained for the hopping t — Q.\Tiq. These Fano-type inter- 
ference objects appear on top of the characteristic spec- 
trum consisting of the Kondo resonance and the broad 
quasiparticle peaks seen at S\ and its Coulomb satellite 
at E\-\-U\. Such spectrum is the result of interference 
effects discussed in the preceding section and the corre- 
lation features. The lower panel of figure |TJ] shows the very 
specific situation £2 = when the Kondo and Fano struc- 
tures coincide with each other. The Fano-type resonance 
seems to play then a dominant role. 

Differential conductance of the Andreev current (Ql for 
the T-shaped double quantum dot system ([I]) is shown 
in figure [71 We notice a clear suppression of the zero- bias 
peak (present for the single quantum dot case t = as 
indicated by the dashed lines) because of a destructive 
influence of the Fano-type interference (the solid lines). 
The subgap Andreev current is thus very sensitive to the 
interplay between the Kondo and Fano effects. For obvi- 
ous reasons their most dramatic competition occurs when 
energy of the side-coupled quantum dot coincides with 
the Kondo resonance, i.e. for £2 = 0. Quantum interfer- 
ence effects destroy then the Kondo peak washing out 
the zero-bias enhancement of the Andreev conductance. 



We have studied a unique nature in which the Fano- 
type quantum interference manifest itself in the energy 
spectrum and differential conductance of the heterojunc- 
tion where a metallic lead is coupled via double quan- 
tum dot to superconducting electrode. In the regime of 
subgap source-drain voltage \eV\ < |A| nonequilibrium 
charge transport is contributed only through the anoma- 
lous Andreev channel when electron from the metallic 
electrode is converted into the Cooper pair (propagating 
further in superconductor) with a simultaneous reflection 
of hole back to the metallic lead. Transmittance of such 
Andreev scattering is a sensitive probe of the proximity 
induced on-dot pairing as well as the quantum interfer- 
ence and correlation effects. 

Since on-dot pairing mixes the particle with hole states 
the interference effects are doubled in a comparison 
to similar junctions without the superconducting elec- 
trode. In particular, for T-shape configuration schemat- 
ically shown in figure [1] we notice that effective spec- 
trum of the interfacial quantum dot develops the res- 
onance and antiresonance, correspondingly at ±£2 (fig- 
ure [5]). These Fano-type structures are present when- 
ever the hopping integral t to the side-attached quantum 
dot (i — 2) is much smaller than the linebroadening Tn 
(whereas Ts merely controls the induced quasiparticle 
splitting). Upon increasing t the Fano-type features dis- 
appear, evolving into the new quasiparticle peaks (figure 
E]) being a consequence of the proximity effect indirectly 
spread onto the side- attached quantum dot. 

Correlation effects play an important role with regard 
to the following aspects: a) the charging effect which 
causes appearance of the Coulomb satellite near £\+U\, 
b) the Kondo singlet state (when the interfacial quan- 
tum dot spin is effectively screened by electrons of the 
metallic electrode leading to formation of a narrow res- 
onance at /j n), and c) eventual suppression the on-dot 
pairing. We have previously shown [161 ] that the Kondo 
effect enhances the zero-bias Andreev conductance as in- 
deed reported experimentally [5] . In the present work we 
indicate that in the double quantum dots the quantum 
interference can (destructively) affect such feature if the 
Fano-type structures appear nearby the Kondo peak. 

A more detailed analysis of the Fano-Kondo interplay 
could be a challenging task in the future studies. For 
this purpose one should resort either to nonperturbative 
techniques, like the numerical renormalization group, or 
to some sophisticated perturbative methods capable to 
interpolate between the limits t — > 0, Tp — > and U — > 0. 
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Appendix: Selfenergy of the noncorrelated DQD 

Using the standard Nambu notation we can express the 
retarded Green's functions of the metallic lead 
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the (unperturbed) side-attached quantum dot 
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and the isotropic superconductor 
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In the equation (IA.3I) we applied the BCS coefficients 
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where £ k =^7+ A 2 . 

For the case of uncorrelated quantum dots (C/j = 0) 
we can determine the selfenergy S°(w) of the interfacial 
quantum dot from the following equation 

E°(«) = J] ^ ffjS(k,«) U^ + t g 2 » t*. (A.4) 

k,/3=7V,S 



Assuming the wide-band limit we introduce the constant 
weighed density of states 



2xJ2\Vw\ 2 H"-tK 



I> for \^ p \<D/2 
elsewhere, 



(A.5) 



where D is the conduction bandwidth. We then easily 
find that 
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because, according to the Kramers-Kronig relation, the 
real part disappears. In the same way we obtain from a 
straightforward algebra that [2( 



X)^Ss(k,w)V k > = 



2« 



(A.7) 



where 



M9(M-A) _ za,e(A-M) fA 



v^ 



VA2 



In the extreme subgap limit \u\ <C A the function (|A.8|) 
approaches j(u>) — > —ico/A and in consequence 



lim S°(w) 
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which proves the equation 
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